Stabilization of Interdomain Interactions in G protein α Subunits as a Determinant of Gαi Subtype Signaling Specificity.

Highly homologous members of the Gαi family, Gαi1-3, have distinct tissue distributions and physiological functions, yet their biochemical and functional properties are very similar. We recently identified PDZ-RhoGEF (PRG) as a novel Gαi1 effector that is poorly activated by Gαi2. In a proteomic proximity labeling screen we observed a strong preference for Gαi1 relative to Gαi2 with respect to engagement of a broad range of potential targets. We investigated the mechanistic basis for this selectivity using PRG as a representative target. Substitution of either the helical domain (HD) from Gαi1 into Gαi2 or substitution of a single amino acid, A230 in Gαi2 with the corresponding D in Gαi1, largely rescues PRG activation and interactions with other potential Gαi targets. Molecular dynamics simulations combined with Bayesian network models revealed that in the GTP bound state, separation at the HD-Ras-like domain (RLD) interface is more pronounced in Gαi2 than Gαi1. Mutation of A230 to D in Gαi2 stabilizes HD-RLD interactions via ionic interactions with R145 in the HD which in turn modify the conformation of Switch III. These data support a model where D229 in Gαi1 interacts with R144 and stabilizes a network of interactions between HD and RLD to promote protein target recognition. The corresponding A230 in Gαi2 is unable to stabilize this network leading to an overall lower efficacy with respect to target interactions. This study reveals distinct mechanistic properties that could underly differential biological and physiological consequences of activation of Gαi1 or Gαi2 by GPCRs.

A novel class of inhibitors that disrupts the stability of integrin heterodimers identified by CRISPR-tiling-instructed genetic screens.

The plasma membrane is enriched for receptors and signaling proteins that are accessible from the extracellular space for pharmacological intervention. Here we conducted a series of CRISPR screens using human cell surface proteome and integrin family libraries in multiple cancer models. Our results identified ITGAV (integrin αV) and its heterodimer partner ITGB5 (integrin ß5) as the essential integrin α/ß pair for cancer cell expansion. High-density CRISPR gene tiling further pinpointed the integral pocket within the ß-propeller domain of ITGAV for integrin αVß5 dimerization. Combined with in silico compound docking, we developed a CRISPR-Tiling-Instructed Computer-Aided (CRISPR-TICA) pipeline for drug discovery and identified Cpd_AV2 as a lead inhibitor targeting the ß-propeller central pocket of ITGAV. Cpd_AV2 treatment led to rapid uncoupling of integrin αVß5 and cellular apoptosis, providing a unique class of therapeutic action that eliminates the integrin signaling via heterodimer dissociation. We also foresee the CRISPR-TICA approach to be an accessible method for future drug discovery studies.

Therapeutic targeting Tudor domains in leukemia via CRISPR-Scan Assisted Drug Discovery.

Epigenetic dysregulation has been reported in multiple cancers including leukemias. Nonetheless, the roles of the epigenetic reader Tudor domains in leukemia progression and therapy remain unexplored. Here, we conducted a Tudor domain-focused CRISPR screen and identified SGF29, a component of SAGA/ATAC acetyltransferase complexes, as a crucial factor for H3K9 acetylation, ribosomal gene expression, and leukemogenesis. To facilitate drug development, we integrated the CRISPR tiling scan with compound docking and molecular dynamics simulation, presenting a generally applicable strategy called CRISPR-Scan Assisted Drug Discovery (CRISPR-SADD). Using this approach, we identified a lead inhibitor that selectively targets SGF29's Tudor domain and demonstrates efficacy against leukemia. Furthermore, we propose that the structural genetics approach used in our study can be widely applied to diverse fields for de novo drug discovery.

Cellular Lipids Regulate the Conformational Ensembles of the Disordered Intracellular Loop 3 in ß2 Adrenergic Receptor.

The structurally disordered intracellular loops (ICLs) of G protein-coupled receptors (GPCRs) play a critical role in G protein coupling. In our previous work, we used a combination of FRET-based and computational methodologies to show that the third intracellular loop (ICL3) modulates the activity and G protein coupling selectivity in GPCRs. In the current study, we have uncovered the role of several lipid components in modulating the conformational ensemble of ICL3 of the ß2-adrenergic receptor (ß2AR). Our findings indicate that phosphatidylinositol 4,5-bisphosphate (PIP2) in the inner leaflet of the membrane bilayer acts as a stabilizing anchor for ICL3, opening the intracellular cavity to facilitate G protein coupling. This interaction between PIP2 and ICL3 causes tilting of ß2AR within the cellular membrane. Notably, this tilting of the receptor is supported by ganglioside GM3 stabilizing the extracellular loops on the outer leaflet of the bilayer, thereby exerting an allosteric effect on the orthosteric ligand binding pocket. Our results underscore the significance of lipids in modulating GPCR activity, proposing an allosteric mechanism that occurs through the receptor's orientation within the membrane.

Bayesian network models identify co-operative GPCR:G protein interactions that contribute to G protein coupling.

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior of interface interactions and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces it is challenging to determine which amino acid pair interactions are cooperative. In this work we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with G proteins. Our results reveal a strong co-dependency in the formation of interface GPCR:G protein contacts. This observation indicates that cooperativity of GPCR:G protein interactions is necessary for the coupling and selectivity of G proteins and is thus critical for receptor function. We have identified subnetworks containing polar and hydrophobic interactions that are common among multiple GPCRs coupling to different G protein subtypes (Gs, Gi and Gq). These common subnetworks along with G protein-specific subnetworks together confer selectivity to the G protein coupling. This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.

Molecular recognition of an aversive odorant by the murine trace amine-associated receptor TAAR7f.

There are two main families of G protein-coupled receptors that detect odours in humans, the odorant receptors (ORs) and the trace amine-associated receptors (TAARs). Their amino acid sequences are distinct, with the TAARs being most similar to the aminergic receptors such as those activated by adrenaline, serotonin and histamine. To elucidate the structural determinants of ligand recognition by TAARs, we have determined the cryo-EM structure of a murine receptor, mTAAR7f, coupled to the heterotrimeric G protein Gs and bound to the odorant N,N-dimethylcyclohexylamine (DMCH) to an overall resolution of 2.9 Å. DMCH is bound in a hydrophobic orthosteric binding site primarily through van der Waals interactions and a strong charge-charge interaction between the tertiary amine of the ligand and an aspartic acid residue. This site is distinct and non-overlapping with the binding site for the odorant propionate in the odorant receptor OR51E2. The structure, in combination with mutagenesis data and molecular dynamics simulations suggests that the activation of the receptor follows a similar pathway to that of the ß-adrenoceptors, with the significant difference that DMCH interacts directly with one of the main activation microswitch residues.

Stabilization of Interdomain Interactions in G protein αi Subunits Determines Gαi Subtype Signaling Specificity.

Highly homologous members of the Gαi family, Gαi1-3, have distinct tissue distributions and physiological functions, yet the functional properties of these proteins with respect to GDP/GTP binding and regulation of adenylate cyclase are very similar. We recently identified PDZ-RhoGEF (PRG) as a novel Gαi1 effector, however, it is poorly activated by Gαi2. Here, in a proteomic proximity labeling screen we observed a strong preference for Gαi1 relative to Gαi2 with respect to engagement of a broad range of potential targets. We investigated the mechanistic basis for this selectivity using PRG as a representative target. Substitution of either the helical domain (HD) from Gαi1 into Gαi2 or substitution of a single amino acid, A230 in Gαi2 to the corresponding D in Gαi1, largely rescues PRG activation and interactions with other Gαi targets. Molecular dynamics simulations combined with Bayesian network models revealed that in the GTP bound state, dynamic separation at the HD-Ras-like domain (RLD) interface is prevalent in Gαi2 relative to Gαi1 and that mutation of A230s4h3.3 to D in Gαi2 stabilizes HD-RLD interactions through formation of an ionic interaction with R145HD.11 in the HD. These interactions in turn modify the conformation of Switch III. These data support a model where D229s4h3.3 in Gαi1 interacts with R144HD.11 stabilizes a network of interactions between HD and RLD to promote protein target recognition. The corresponding A230 in Gαi2 is unable to form the "ionic lock" to stabilize this network leading to an overall lower efficacy with respect to target interactions. This study reveals distinct mechanistic properties that could underly differential biological and physiological consequences of activation of Gαi1 or Gαi2 by GPCRs.

Dynamic spatiotemporal determinants modulate GPCR:G protein coupling selectivity and promiscuity.

Recent studies have shown that G protein coupled receptors (GPCRs) show selective and promiscuous coupling to different Gα protein subfamilies and yet the mechanisms of the range of coupling preferences remain unclear. Here, we use Molecular Dynamics (MD) simulations on ten GPCR:G protein complexes and show that the location (spatial) and duration (temporal) of intermolecular contacts at the GPCR:Gα protein interface play a critical role in how GPCRs selectively interact with G proteins. We identify that some GPCR:G protein interface contacts are common across Gα subfamilies and others specific to Gα subfamilies. Using large scale data analysis techniques on the MD simulation snapshots we derive a spatio-temporal code for contacts that confer G protein selective coupling and validated these contacts using G protein activation BRET assays. Our results demonstrate that promiscuous GPCRs show persistent sampling of the common contacts more than G protein specific contacts. These findings suggest that GPCRs maintain contact with G proteins through a common central interface, while the selectivity comes from G protein specific contacts at the periphery of the interface.

Cytochrome c oxidase inhibition by calcium at physiological ionic composition of the medium: Implications for physiological significance of the effect.

Cytochrome c oxidase (CcO) from mammalian mitochondria binds Ca2+ and Na+ in a special cation binding site. Binding of Ca2+ brings about partial inhibition of the enzyme while Na+ competes with Ca2+ for the binding site and protects the enzyme from the inhibition [Vygodina, T., Kirichenko, A. and Konstantinov, A.A. (2013). Direct Regulation of Cytochrome c oxidase by Calcium Ions. PLoS One 8(9): e74436]. In the original studies, the inhibition was found to depend significantly on the ionic composition of the buffer. Here we describe inhibition of CcO by Ca2+ in media containing the main ionic components of cytoplasm (150mM KCl, 12mM NaCl and 1mM MgCl2). Under these conditions, Ca2+ inhibits CcO with effective Ki of 20-26µM, that is an order of magnitude higher than determined earlier in the absence of Na+. At physiological value of ionic strength, the inhibition can be observed at any turnover number of CcO, rather than only at low TN (<10s-1) as found previously. The inhibition requires partially oxidized state of cytochrome c and is favored by high ionic strength with a sharp transition at 0.1-0.2M. The high Ki=20-26µM found for CcO inhibition by calcium matches closely the known value of "Km" for Ca2+-induced activation of the mitochondrial calcium uniporter. The inhibition of CcO by Ca2+ is proposed to modulate mitochondrial Ca2+-uptake via the mitochondrial calcium uniporter, promote permeability transition pore opening and induce reduction of Mia40 in the mitochondrial intermembrane space.